extern int drotg_(doublereal *da, doublereal *db, doublereal *c__, doublereal *s);
